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Abstract 

The intensity of electron spin resonance (ESR) of the nanoscale molecular magnet V15 is stud- 
ied. We calculate the temperature dependence of the intensity at temperatures from high to 
low. In particular, we find that the low-temperature ESR intensity is significantly affected by the 
Dzyaloshinsky-Moriya interaction. 
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I. INTRODUCTION 



The Vi5 molecule has been one of promising nanometer-scale molecular magnets since 
it was first synthesized.-"- It is the complex of formula Kg 
Vi5, fifteen vanadium ions of spin 1/2 form almost a sphere. Three spins in the middle are 
sandwiched by the upper and lower hexagons. 

Different experiments on the magnetization process have shown that the magnetization 
changed adiabatically in a fast sweeping field, and a magnetic plateau appeared in a slow 
sweeping field due to thermal bath attached to the molecule.--^ The latter phenomenon, 
which is called the phonon bottleneck effect, is theoretically analyzed from a general point 
of view of the magnetic Foehn effect.- This smooth change of the magnetization at H = OT 
implies the existence of an avoided-level-crossing energy structure. The structure of avoided 
level crossing has been studied. In a model of the triangle Heisenberg antiferromagnet with 
three spins, at if = 0, two sets of 5* = 1/2 doublets overlap, and 4 states degenerate. The 
degeneracy is resolved into two sets of Kramers doublets by perturbation such as anisotropy, 
Dzyaloshinsky-Moriya interaction (DMI),— "— and also the hyperfine interaction.— Depend- 
ing on the type of perturbation, there appear many kinds of energy structure. At the 
crossing point of the two states of m = 1/2 of the doublet states and one of the quartet 
state (m = 3/2), H = He, a kind of avoided level crossing is formed. It has been pointed 
out that adiabatic change causes the change of magnetization from m — l/2tom = 1 be- 
cause one state adiabatically changes to the state with m = 3/2.— In V15, the equilibrium 
magnetization curve shows smooth change at zero field from 1/2 to 1/2 and at 2.8T (= He) 
from 1/2 to 3/2. V15 can be described by a triangle model but details of the DMI in V15 
are not yet fully understood. 

In this paper, first we numerically calculate the temperature dependence of the ESR 
intensity of V15 using a new numerical method (the double Chebyshev polynomial method) 
of calculating the Kubo formula. We find that the model Hamiltonian for V15 including the 
DMI successfully reproduces the experimental temperature dependence of the ESR intensity. 
Second we investigate the ESR at very low temperatures. We find that the intensity ratio 
(the intensity of V15 divided by that of a spin 1/2) is affected by the DMI at small fields. 
We propose that experimental observation of the intensity ratio enables us to estimate the 
DMI in V15. Finally, we analyze the ESR at low temperatures using a triangle model whose 
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energy levels model the low-lying levels of V15. 
II. MODEL AND FORMULATION 
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FIG. 1: Schematic picture of V15. 



Figure [T] shows the structure of vanadium ions in V15. We consider the following spin 
Hamiltonian for V15.— 



The first term on the right-hand side of Eq. ([T]) describes the Heisenberg interaction. Coef- 
ficients Jij take three values J, Ji, and J2 > I -^2 1 > depending on the bonds on the 
upper and lower hexagons. Three spins between two hexagons interact with the hexagons 
by Ji and J2. There is no interaction among these three spinsS'. We set J = — 800K, 
J2 = -350K, and Ji = -225K.-^ The second term describes the DMI. We assume the 
existence of DM vectors at the bonds of J. In the third term, H denotes the static 

magnetic field applied to the molecule. We will ignore other effects such as dipolar fields, 
hyperfine interactions, and the crystal field, which are considered to be negligibly small. 
Indeed, the dipolar and hyperfine fields are estimated as 1 mK and 50 mK, respectively^. 
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FIG. 2: Heisenberg interactions between the spins. 



Figure E] explains the interactions between spins. If we assume the symmetry of V15, 
the lower hexagon differs from the upper hexagon by rotation 7r/6 and we have only one free 
DM vector, say Di^2 (Fig.E]). We take the vector Di^2 to be D^^ = Df ^ = ^^,2 = 40K. Let 
us define 
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We obtain the other DM vectors on the upper hexagon by rotating -Di,2 by 27r/3 and Att/3: 

£>3,4 = R (^] D,^2, D,^, = R C^) £)i,2, (3) 



i.e., 4 = 14.641K, = -54.641K, D^^ = 40K, g = -54.641K, e = 14.641K, and 
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40K. The DM vectors on the lower hexagon are obtained as 
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Let us calculate the ESR intensity of V15 using the Hamiltonian ([T]). By the Kubo 
formulai^ii^^, the imaginary part of the susceptibility ^("{(jJ^T) on the transverse field H 
parallel to the x-axis is given by the Fourier transform of the spin-spin correlation function; 

T) = (1 - e"^^) Re {M^M^{t))e-'"^'dt, (5) 

where u is the frequency of the radiation field, T is temperature, and /3 = 1/T. The 
absorption is given by 

/(^;T) = ^x"(^;T), (6) 

where i^R is the amplitude of the radiation field. We obtain the ESR intensity /(T) by 
integrating J(a;,T) with respect to uj] 

poo 

I{T)= I{u,T)du. (7) 
Jo 

In the present paper, we obtain I{T) without diagonalization. Our method can readily 
be applied to other nanomagnets. In Appendix \^ we describe the numerical method of 
calculating the Kubo formula for huge Hamiltonian quantum many-body systems.-^S This 
method, which we call the double Chebyshev expansion method (DCEM), realizes 0{N) 
calculation in both speed and memory. In the method, we evaluate the summation over 
all the states in the expression of the correlation function by using the average over a few 
number of random initial states. Furthermore, we calculate the exponential operators e~^^ 
and e~^^* by expanding them in the Chebyshev polynomial. The DCEM takes all the states 
in the Hilbert space into account and enables us to obtain the ESR at any temperature. 



5 



FIG. 3: The D3 symmetry of V15. 



III. TEMPERATURE DEPENDENCE OF THE ESR INTENSITY 

We study the temperature dependence of the ESR intensity of V15. We apply the magnetic 
field parallel to the c-axis of the molecule (z-axis): H = (0, 0, H). Figure S] shows intensities 
at i^" = 2T from IK to lOOOOK. The intensity obtained by the DCEM (solid circles) is 
plotted together with experimental data by Ajiro, et al.^- (solid squares). The intensity by 
the DCEM in Fig. H] agrees with the experimental data. In addition, the dashed line denotes 
the intensity 

h{T) = Htanh^ (8) 

of an isolated spin 1/2. The short-dashed line, dotted line, and dash-dotted line denote 
2/1 (T), 3/1 (T), and 15/i(T) for non-interacting 2, 3, and 15 spins, respectively. 
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FIG. 4: The temperature dependences of the intensity for H\\z and H = 2T and the experimental 
intensity. 

When the temperature is sufficiently higher than the interactions, spins in V15 are almost 
isolated. Therefore, the intensity is expected to meet the dash-dotted line at very high 
temperatures. In Fig. HI the intensity starts to deviate from the dotted line near 200K. 
As the temperature decreases, the effective number of spins changes from 15 to 3, and the 
intensity stays on the dotted line at temperatures between lOK and lOOK. 
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FIG. 5: The temperature dependence of the intensity for H\\z and H = 4T. 



Figure |5] is the same as Fig. H] except H = 4T. At temperatures above lOK, the intensity 



in Fig. [5] behaves almost the same as that in Fig. HI However, the temperature dependence 
of the effective number is different from that in Fig. H] below lOK. In Fig. |5l the effective 
number of spins changes from 15 to 3 at high temperatures, and stays on the dotted line for 
3/i(T) as temperature decreases. A similar behavior of the low-temperature intensity has 
also been predicted.— 



IV. THE LOW-TEMPERATURE INTENSITY RATIO 



A. I{T) at low temperatures 

Although the DMI affects I{T) only mildly at T > IK, it significantly changes /(T) at 
low temperatures. At low temperatures, only transitions among low-lying energy levels near 
the ground state are responsible for the ESR absorption. This fact enables us to obtain 
the intensity by investigating the transitions among the lowest eight levels. Thus, we can 
calculate intensity with the subspace iteration method (SIM).— "^^i^^ The response x"{^yT) 
is obtained by direct diagonalization in the small reduced space. In Fig. [6l we compare the 
result of Fig. H] and the data obtained by the SIM. We find a good agreement below lOOK. 
Above this temperature, the lowest eight states are not enough to represent the system. 
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FIG. 6: The temperature dependence of the intensity for H\\z and H = 2T. 

In Fig. [Tl^a), we show the lowest eight energy levels of V15 as a function of the magnetic 
field H applied in the z-direction. Due to the DMI, we have an avoided level crossing near 
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3T as shown in Fig. Wijo)- Without the DMI, the gap closes and the two levels just cross 
each other. 




FIG. 7: (a) The lowest eight levels of V15 as a function of H\\z. (b) Magnified figure near 3T. 

We define the intensity ratio R(T) as 

RiT) = /(T)//i(T), (9) 

where Ii{T) is the intensity of a single spin (Eq. ([8])). In the previous section, we have found 
R{T) ~ 3 for if = 4T and R{T) is slightly larger than 1 for H = 2T in the low temperature 
limit. Figure m shows R{T) in a low temperature region. In Fig. [HI the circles and triangles 
show R{T) at H = 57.8GHz (~ 2T) and at H = 108GHz (~ 4T), respectively. The dashed 
lines in Fig. [8] show the ratios without the DMI at H = 2T and at H = 4T for comparison. 
First, we notice that He can be experimentally estimated with the low-temperature ESR by 
measuring the field at which the destination of the intensity changes. Second, we see that 
R{T = OK) at ii = 2T deviates from 1 due to the DMI, while R{T = OK) at H = AT stays 
very close to the dotted line. Thus it would be possible to experimentally determine the 
DMI in Vi5 by observing the intensity at low temperatures. 
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FIG. 8: The intensity ratios of V15 as a function of the temperature. The intensity ratios without 
the DMI are also shown for comparison. 

Furthermore, we study intensity ratios at various fields near the avoided level crossing 
point (hereafter we refer to this field as He) in the presence (Fig.[9](a)) and absence (Fig.[9]^b)) 
of the DML The ratio goes to 3 as T —t- for a field higher than H^. In contrast, it goes 
near 1 for a lower field. The derivative of the line for H = 2T in Fig. [9](a) is negative at OK, 
while that for if = 3T is positive. Hence the line for H = 2T has a dip at T ~ 0.5K. In the 
absence of the DMI, the line for H = 2T does not have a dip. 
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FIG. 9: (a) Intensity ratios of V15 as a function of temperature for various fields, (b) Intensity 
ratios without the DMI. 



10 



B. Triangle Model 



As shown in Figs. [8] and [9]^a), the intensity ratio for lower fields deviates from 1 in the 
low-temperature limit. Here, we study the mechanism of this deviation analytically in a 
triangle model of three spins which well describes the low-lying energy structure of V15. 

The triangle modeU^ is given by the Hamiltonian Eq. ([1]) with Jjj = J < 0, H = {0,0, H), 
H>0, 

~ = Dy, 

Dh = i-D. + V^Dy)/2, Dl^ = {-VSD, - Dy)/2, (10) 
D^i = i-D^ - VSDy)/2, Dl = {Vsa, - Dy)/2, 

and 

Dl, = D^,, = Dl, = D^. (11) 

Note that the elements of the DM vectors are chosen so that the model satisfies the C3 
symmetry. If we put J = — 2.5K and Dj. = Dy = = 0.25K(= D), energy levels of the 
triangular model well reproduce the lowest eight levels of V15 as shown in Fig. [101 In the 
absence of the DMI, the ground state magnetization changes from 1/2 to 3/2 at 

H, = ^\J\ ^2.8T. (12) 




FIG. 10: Energy levels of the triangle model as a function of the field H. 

The intensity ratio RaiT) of this triangle model is obtained as 
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r(m) 

m 

J2 {En - Em) (e-^^™ - e-^^") |2 {Wm\ M,. \Wr,)f , (13) 

where {Em} and {|u^m)} are the eigenvalues and eigenvectors of the triangle-model Hamil- 
tonian, respectively. In Appendix [Bl RsiT) at low temperatures around He is calculated. 

Let us explore the nonmonotonic temperature dependence of the intensities shown in 
Fig. [9] by using the triangle model with different values of D^, Dy, and D^. 

First, we consider how the DMI causes the dip at T ~ 0.5K for H = 4T in Fig. [9]^a). In 
the absence of the DMI (Fig.^b)), by setting D^ = Dy = D, = 0, we obtain dR^{T)/dT > 
for H < He and di?3(r)/dr < for H^ < H (see Eq. (!BT5|) in Appendix E]) . To investigate 
the intensities in Fig. IHt^a), we set = Dy = 0, ^ (uniaxial). Note that the structure 
of the energy levels with Dr^ = Dy = D^ ^ is the same as that in the uniaxial case only 
except that the gap closes in the uniaxial case. We find that dR3{T)/dT < for H ^ H^ 
and di?3(T)/dT > for if ~ ifc (see Eq. (1B17I) in Appendix iBl). which implies a dip at a 
low field. 

Next, let us consider the effect of a weak DMI. We set D^ = Dy = D^ = D ^ 0. By 
ignoring smaller terms than 0{D), we have (Appendix [B|) 

I 3 _ (^2 - I) e-l^-^-^/2|/^ (He <H), 

where we defined 

A = V3D. (15) 

Taking the limit T — )■ 0, we have 

f 1 + # (H <He) 
R,{T) ^ ^ ^ (16) 

[3 {He< H). 

This term A/if describes the deviation of R{T) from 1, and the deviation is proportional 
to D. 



V. SUMMARY 



In this paper, we studied the ESR of the nanomagnet V15. We investigated the temper- 
ature dependence of the intensity on the DMI for different values of the static field H. The 
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DMI significantly affects the low-temperature intensity of V15. We found that the intensity 
at H < He has a dip as a function of temperature due to the DMI. We analyzed the dip us- 
ing the three-spin model and obtained analytical expressions of the intensity. Experimental 
measurement of the dip may be used to estimate the DMI of V15. 
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Appendix A: DCEM 

The DCEM is an extension of the Boltzmann-weighted time-dependent method 
(BWTDM).-21>28 The DCEM differs from the BWTDM by the third step below. The proce- 
dure of the DCEM is divided by the following five steps. 

At the first step, we prepare a random vector |<I>).i^i^ For a given basis \n) of the Hilbert 
space, this random vector is given by |$) = Y,n=i l'^) Here, the dimension of the Hilbert 
space is N and the random numbers {On} take values from — vr to vr. 

At the second step, we obtain the Boltzmann-weighted random vector |$Boitz) = 
|(|,^ expanding it with the Chebyshev polynomial; 

e-/5«/2 = Jo(-/3AA/2)ro(?^sc) 

+ 2j2lk{-f3AX/2)Tk{'Hsc), (Al) 

k 

where Ik{x) is the modified Bessel function and Tfc('Hsc) is the Chebyshev polynomial, which 
satisfies Tk{V.sc) = 2'HscTk-i{'Hsc) - Tk-2{'Hsc), To('Hsc) = 1, and Ti('Hsc) = Hsc- We note 
that the eigenvalues of "Hsc (= "H/AA) are confined between —1 to 1. By this procedure, we 
obtain vectors e~^^/^ |$) without diagonalization. 

At the third step, we obtain |$Boitz;i^) = e~^^* |$Boitz) and \^M^',t) = e~'^* \^m^), where 
= M"" |$Boitz) and M"" = Sj. In the BWTDM, the time evolution is performed by 
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the leap frog method, which evolves a state as 

10; t + At) = -2inAt 10; t) + \(j);t- At) . (A2) 

Note that the condition E^aaxAt <^ 1 should be satisfied, where -Bmax is the largest eigenvalue 
of the Hamiltonian. In the DCEM, we make use of the Chebyshev polynomial expansion in 
order to obtain the time evolution of the vector; 

\(f);t + r) - e-^"^^«-|0;i) 

= MTAX)To{nsc)\<P;t) 

+ 2 ^(-i) Vfe(rAA)Tfc(?^,e) 10; t) , (A3) 

k 

where Jk{x) is the Bessel function. Note that the time step r is not necessarily small. In 
the ESR experiment for V15, the magnetic field H (~ IK) is usually much smaller than the 
strongest couphng \J\ (~ lO^K). Hence the frequency of precession of the spins is small. This 

means that we need to evolve state vectors for a long time but do not need fine resolution 
of the time step. This is why the DCEM is more efficient than the BWTDM for the ESR of 
V15. 

At the fourth step, we calculate the correlation function 



(M^M^it)) - 

[($M^;t|M-|$Boitz;t)],v ^^^^ 



Jav 



[(^Boltzl^Boltz)] 

where the traces are replaced by the ensemble averages ([-jav) with respect to the random 
vectors; for any operator X, TrX is calculated as 

n m,n 

~ TrX. (A5) 

Finally, x"(a;;T) is obtained by the Fourier transform of {M^M^{t)). 

x"(cu;T) = (1 - e-'^") Re (M=^M^(t))e-^'"*dt 

1 - e"'^'") Re y (M^M^(t))e-^'"*e-^ * /Mi. (A6) 

Here we introduced the Gaussian filter with variance l/r;^. This 77 determines the frequency 
resolution. The upper limit of the integral Tmax satisfies Tmax ~ I/77 in order to avoid 
the Gibbs oscillation. Also 77 should satisfy the conditions that < rj <^ 1, rj <^ H, and 
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Appendix B: Details of the Triangle Model 

We obtain the block-diagonalized form of the Hamiltonian by using the following basis 
vectors 

, , -1 



\V2) = 

1^4) = 
1^5) = 



2^/3 
-1 

-1 

-1 

\m) 



(1 + iV3) im) + (1 - iV3) \m) - 2 \m) 
(1 - ix/3) im) + (1 + i\/3) \m) - 2 i;n) 
-2 im) + (1 + iV3) im) + (1 - iv^) i;;t) 
-2 im) + (1 - iVs) \m) + (1 + iV3) lut) 



^ [im) + int) + i;tt)] , 



l^;^) = -j= [\m) + im) + m)] , and 



(Bl) 



Now we have 



n\v4) 

n\v2) 
n\vr) 

nivs) 



= -l{3J + 6H)\v,) + lD+\v,), 

= |D_ \V5) + \{3J-2H- 2x/3D,) \v,) , 

= |(3J + 2// + 2V3D,) \v,) + \ve) , 

= fD^\v,)-\i3J + 2H)\v,), 

^ \{3J-2H + 2V3D,) \V2) - \v,) , 

= -^D^\v2) + l{-3J + 2H)\vj), 

= i(-3J + 6f/)K) + |D_|^3), and 
= |L>+ |^;8) + 1{3J + 2H- 2y/3D,) \v^) , 



where D± — iDy. 



(B2) 



(B3) 



(B4) 



(B5) 



a. Uniaxial Dzyaloshinsky-Moriya interaction 



In the presence of the uniaxial DMI [D^ = Dy = 0), the triangle-model Hamiltonian is 
diagonalized by the vectors {l^j)}, i.e., {I^j)} form the eigenvectors of the Hamiltonian. The 
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correspondent eigenvalues {Ej} are given by 



Eg — 


3 3 
— 7+ _ff 
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^7 — 


4-^ + 2^' 




Eq — 


-'-J-\ 

A ' 




E, = 






E, = 


"4: 




Es — 


4 2 
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E2 = 


4 2 


—^Dzi and 


El = 


4 2 


2 ^ 



Besides, nonzero matrix elements are obtained as 



(^1 


\Vi) 




\v:,) 




\ve) 




\V7) 


{V7 


\vs) 



1 

~2' 
1 

~2' 

y5 

2 ' 

1, and 
2 ' 



where 



= + SI + SI. 

Note that iVP' is hcrmitian, i.e., {vi\ \vj) = {{vj\ \vi))*- 
The magnitude relation of Ei, E2, and E^ depends on H: 

' Ei<E2<E, {H<H,-f), 
Ei<E^<E2 {H,-^<H <H, + ^), and 
E^<Ei<E2 {H, + ^<H), 

where we define 
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In the absence of the DMI {D^ — Dy — Dz — 0), the two- fold degenerate ground state 
(\vi) , \v2)) of Mz = 1/2 and the state (Iv^)) of = 3/2 intersect at if = H^. In the 
uniaxial DMI, the states 1^2) and {v^) intersect at H — He — A^, and the states \vi) and l^s) 
intersect at H = He + A^, respectively. 

Note that, near H^, we calculate the intensity ratio by taking into account only three 
low-lying states \vi), \v2), and 1^5) at very low temperatures, and we have 



-l3Er, 



+ e' 



^r(m) ~ r(l) -I- r(2) -I- r(5) 

m 

~ (// + A^)e-'^^i + He-^^^ + 3He~^^'. 
Hence RsiT) is approximated as 



tanh (f^) (e-z'^i + e-/^^^ + q-^e,^^ 



3 + g-/3(i/-i/c+A,/2) 



1 + g-/3(H-Hc-A./2) 



1 + g-/3(i/-i/c+A^/2) _^ g-/3(H-//c-A^/2) 



In the absence of the DMI (A^ = 0), we have 

~ 1 + 2e-/3(^^-^c) ■ 
Therefore, at ultra-cold temperatures (T <^ |if — i^d), we have 

1 + e-\H-H.\/T < 



and 



i?3(r) ~ 



dRsiT) 



3 _ 4e-l^-^^c|/T 



In the case of finite A^, at ultra-cold temperatures {T <^ A^ for H < H^ 
T^\H-He- AJ2\ for H^ - AJ2 < H), we have 



(Bll) 



(B12) 
(B13) 

(B14) 

(B15) 
A^/2 and 



1 I A, _ 2A, „-A,/T 
^ ^ H H ^ 



{H<He- A,/2) 
1 + ^ + (2 - ^) Q~W~H.-^./2\/T ^ff^ _ <H<He + AJ2) 

3 _ (2 - ^) Q-\H-H.-^./2\/T ^ ^^/2 < H), 



(B16) 
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and 



dRs{T) 
dT 



— < 



_i_ 

T2 



9 A^ 



\H-Hr- 



{H<H,- A,/2) 
g_|H-ff,-A./2|/T _ <H<H, + A,/2) 



(B17) 

In the uniaxial case, the structure of the energy levels is almost the same as that of V15 
except that the gap at the avoided level crossing point closes. The equation explains the 
behavior of R{T) in the simulation of V15: R'{T) < at low fields, R'{T) > at fields close 
to the avoided level crossing, and R'{T) < at high fields. 



b. Weak Dzyaloshinsky-Moriya interaction 

Let us consider the full DMI case {D^ = Dy = = D). We assume D > is small in 
the sense that D <^2\H — Hc\ is satisfied. 

By ignoring smaller terms than 0{D), we have the following eigenvectors {|%)}. 

\ui) = -V3e-^-/^e ^5) + l^^i) 

\U2) = bs) + VSe'^/^^ \vi) 

\us) = -e-'^/^rj \V2) + \vr) 

\u,) = -VSe^^/^r? |^;8) + |^;3) 

\U5) = \v2)+e'''/%\v7) 

\U7) = e'^/^^\v,) + \ve) 

\us) = \vs) + V3e-'''/^r]\v3), 

(B18) 

where 

e = — , 77 = — . (B19) 

The corresponding eigenvalues {Ej\ are given by 

18 







1^ 
— — ij 

2 




P — 


4"^ 


3^ 
2 




P — 




2 + 


A 


TP 

-C/3 — 


4^+ 


1^ 
— jj — 
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El = 
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A 

Y 



(B20) 



Prom the above eigenvectors {]%)}, nonzero matrix elements are obtained as 



(til 




1 ^4/ 


1 

~ ~2 


{Ui 








{U2 




1^3) 


1 

~ ~2 


{U2 




\ua) 


= 


{U2 




\uq) 


= e-'^/^77 


{U2 






= 


{U3 




1^7) 




{U4 






= -v^e'^^/^e 


(m4 




1^7) 


= -e^"/^e 


{U5 






_ x/3 
2 


{U6 






= 1 


{U7 






_ x/3 
2 ■ 



Let us consider the uhra-cold hmit (/3 — )■ 00), where aU transitions can be ignored except 
for the transitions from the ground state. We consider two cases: H < H^, where l^i) is the 
ground state, and < H, where 1^5) is the ground state. We have 

m 

J2 r(m) ~ r(l) + r(2) ~ (// + A) e'^^^ (H - A) e"^^^ 

m 

(B22) 
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for H < Hr and 



r(m) ~ r(5) + r(l) ~ SHe'^^' + {H + A) e'^^' 

m 

for He < H. Thus we obtain Eq. (fTil) . 



(B23) 
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